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The  growth  of  high-quality  organic  semiconductors  on  sup¬ 
porting  substrates  with  minimal  disorder  is  the  key  require¬ 
ment  for  realizing  high-performance  organic  electronics  J1’2] 
However,  the  absence  of  an  epitaxial  relation  and  the  presence 
of  disorder  make  it  particularly  challenging  to  grow  high-quality 
organic  crystals  having  few  defects  and  grain  boundaries  on 
conventional  insulating  substrates. Here,  we  report  epitaxial 
growth  of  highly  ordered,  crystalline  organic  films  on  van  der 
Waals  (vdW)  substrates  and  demonstrate  high  carrier  mobility 
comparable  to  those  of  single-crystal  counterparts.^  High- 
quality  rubrene  films  with  large  single-crystalline  domains 
were  grown  on  hexagonal  boron  nitride  (h-BN)  layers  with 
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an  in-plane  crystallographic  registry.  Furthermore,  using  gra¬ 
phene  as  the  vdW  electrical  contacts,  high  field-effect  mobility 
(11.5  cm2  V-1  s-1)  was  achieved,  attributable  to  an  atomically 
sharp  interface  with  low  charge  trap  density  between  rubrene 
and  h-BN.  The  ability  to  grow  high-quality  organic/inorganic 
epitaxial  vdW  hetero structures,  combined  with  recent  progress 
on  large-area  growth  of  layered  materials, provides  new 
opportunities  for  the  scalable  fabrication  of  high-performance 
organic  electronic  devices. 

Controlling  the  interface  between  organic  films  and  inor¬ 
ganic  substrates  is  of  crucial  importance  for  organic  electronic 
applications,  because  this  interface  governs  the  initial  nucleation 
during  the  growth  of  molecular  films  as  well  as  the  field-induced 
charge-carrier  transport^2’8-101  In  particular,  charge  transport  is 
affected  greatly  by  interfacial  disorder,  such  as  surface  roughness 
and  charge  traps,  as  well  as  structural  disorder  of  the  molecular 
materials  at  the  interface  J9,11-14]  The  highest  carrier  mobility 
and  intrinsic  band-lilce  transport  have  been  achieved  mostly 
with  organic  single  crystals  grown  by  a  physical  vapor  trans¬ 
port  method  and  assembled  onto  surface-functionalized  device 
structures  or  suspended  over  an  “air-gap” J15,16]  Although  this 
approach  provides  a  platform  for  fundamental  transport  studies, 
the  complicated  device  geometry  and  fabrication  processes  pre¬ 
vent  its  practical  implementation.  Therefore,  it  is  desirable  to 
grow  high-quality  crystalline  films  directly  on  functional  sub¬ 
strates. [17_2°]  However,  device  performance  comparable  with  that 
of  free-standing  single-crystal  devices  has  not  been  achieved  due 
to  the  difficulty  in  forming  organic  films  with  long-range  struc¬ 
tural  ordering  on  typical  substrates  such  as  silicon  oxide,  glass 
and  plastic,  which  possess  a  high  degree  of  disorder. 

Unlike  amorphous  substrates,  two-dimensional  (2D)  dielec¬ 
tric  materials  such  as  h-BN  and  other  vdW  materials  such  as 
metal  oxides  and  metal  chalcogenides  offer  a  new  type  of  sub¬ 
strate  for  high-quality  epitaxial  growth  of  organic  semiconductor 
films.  In  general,  the  2D  crystal  lattice  of  layered  materials 
should  allow  organic  molecules  to  assemble  into  well-ordered 
structures  via  vdW  interactions  while  minimizing  random  het¬ 
erogeneous  nucleation  processes. In  addition,  a  vdW  surface 
is  free  of  dangling  bonds  and  charge  traps  and  thus  enables  the 
formation  of  the  atomically  well-defined  interface  with  organic 
materials.^21!  Reducing  the  density  of  the  charge  traps  and/or 
scattering  centers  is  desirable  for  high-performance  devices. 

In  this  study,  we  demonstrate  the  growth  of  high-quality 
rubrene  crystals  on  exfoliated  h-BN  crystals.  Figure  la  contains  a 
schematic  of  a  rubrene/h-BN  heterostructure  after  crystal  growth. 
It  is  important  to  note  that  the  nonplanar  molecular  structure  of 
rubrene  favours  assembly  on  the  h-BN  surface  with  an  “edge-on” 
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Figure  1.  Growth  of  organic/2D  layered  vdW  heterostructures.  Schematic  illustration  of:  a)  the  organic/2D  layered  vdW  heterostructure  composed  of 
rubrene  and  h-BN  and  b)  the  vapor-phase  transport  method  for  the  growth.  Optical  images  of:  c)  the  h-BN  crystal  exfoliated  on  a  Si02/Si  substrate  and 
d)  the  rubrene  film  grown  on  h-BN.  e)  AFM  image  of  the  film  surface  and  the  height  profile  along  the  line.  It  shows  a  highly  uniform  terrace  structure 
with  the  step  height  of  1 4  ±  1  A. 


orientation  rather  than  “face-on”  alignment.  This  arrangement 
offers  an  advantage  for  lateral  charge  transport  parallel  to  the  semi¬ 
conductor/dielectric  interface  in  a  field-effect  device  geometry.^22’23! 
Figure  lb  shows  a  schematic  of  the  vapor-phase  transport  method 
used  for  the  growth  of  rubrene  crystals  (see  Experimental  Section 
for  details).  To  prepare  the  substrate  for  growth,  we  created  thin 
layers  of  h-BN  crystals  that  were  mechanically  exfoliated  on  Si02/Si 
substrates  as  shown  in  Figure  lc.  Typically,  single  h-BN  crystallites 
have  a  thickness  of  10-60  nm  and  a  lateral  dimension  larger  than 
10  pm.  We  chose  high-quality  crystallites  that  have  extremely  flat 
surfaces  (roughness  of  less  than  1  nm)!24'  in  order  to  provide  atom¬ 
ically  homogeneous  substrate  for  the  subsequent  organic  crystal 
growth.  After  cleaning  the  h-BN  surface  via  thermal  annealing, 
homogeneous  rubrene  films  with  thicknesses  of  500-1000  nm 
were  preferentially  grown  on  the  h-BN  surface  (Figure  Id)  along 
with  microcrystallites  occasionally  formed  on  the  Si02  substrate  at 
the  optimal  growth  condition  (see  Supporting  Information,  Figure 
SI).  The  atomic  force  microscopy  (AFM)  image  of  the  rubrene 
on  h-BN  shows  atomically  flat,  highly  uniform  terraces  with  an 
average  step  height  of  14  ±  1  A  (Figure  le).  The  step  height  cor¬ 
responds  to  the  molecular  length  of  the  c-axis  of  rubrene. 

Figure  2a  shows  a  transmission  electron  microscopy  (TEM) 
image  and  selected  area  electron  diffraction  (SAED)  pattern  of 
the  rubrene  film  grown  on  an  h-BN  crystal.  For  this  measure¬ 
ment,  we  place  thin  h-BN  crystals  on  a  holey  TEM  grid,  and 
then  grow  the  rubrene  crystals.  The  rectangular  diffraction 
patterns  from  rubrene  (hkO)  planes  are  clearly  seen  on  top  of 
the  hexagonal  diffraction  patterns  of  h-BN  {lOlo}  planes.  The 
superposition  of  the  two  diffraction  patterns  without  concentric 
rings  clearly  indicates  that  both  h-BN  and  rubrene  layers  are 
crystalline  and  no  amorphous  layer  is  formed  at  the  interface 
between  the  two  crystals.  Furthermore,  by  comparing  the  2 -fold 
diffraction  pattern  of  rubrene  to  the  6-fold  symmetric  pattern 
of  h-BN,  we  find  that  the  two  single  crystals  are  rotated  by  an 
angle  of  ca.  4°  between  their  a-axes. 


Since  SAED  can  only  provide  local  structural  information,  we 
used  grazing-angle  X-ray  diffraction  (GIXD)  with  a  synchrotron 
light  source  (Figure  2b).  Due  to  the  highly  crystalline  nature  of 
the  films,  we  observe  sharp,  intense,  and  well-defined  Bragg  dif¬ 
fraction  peaks.  Only  a  small  fraction  of  the  total  diffraction  peaks 
were  observed  at  a  given  0  angle.  The  diffraction  peaks  appear 
and  disappear  over  a  narrow  angular  0  range  (see  Supporting 
Information,  Figure  S2).  We  thus  obtain  a  complete  reciprocal 
(g-space)  diffraction  pattern  by  measuring  while  the  sample  is 
rotated  360°  around  the  film  normal  axis  as  shown  in  Figure  2c. 
The  experimental  data  are  in  a  good  agreement  with  the  com¬ 
puted  diffraction  patterns  (indicated  by  green  circles)  from  an 
orthorhombic  unit  cell  with  the  lattice  parameters  (a  =  7.19,  b  = 
14.43,  c  =  26.83  A).[25]  This  reveals  c-axis-oriented  growth  of  an 
orthorhombic  rubrene  crystal  on  the  h-BN  substrate.  The  lat¬ 
tice  parameters  obtained  from  GIXD  are  close  to  the  values  of 
a  free-standing  single  crystal,  indicating  that  the  intermolecular 
interactions  in  our  rubrene  film  are  stronger  than  the  interac¬ 
tion  between  the  rubrene  molecules  and  the  h-BN  substrate.[21] 

The  agreement  between  experimental  and  computational  dif¬ 
fraction  patterns  enables  a  complete  peak  indexing  and  deter¬ 
mination  of  in-plane  crystallographic  orientation  of  the  rubrene 
film  relative  to  the  h-BN  substrate.  Figure  2d  plots  the  intensity 
as  a  function  of  0  for  a  given  diffraction  peak.  Interestingly,  we 
find  that  rubrene  crystals  grown  on  h-BN  have  6-fold  symmetry 
even  though  the  orthorhombic  crystal  structure  of  rubrene  only 
has  2-fold  symmetry.  Considering  that  there  is  a  possibility 
of  multiple  domains  of  rubrene  on  h-BN,  the  observed  6-fold 
GIXD  suggests  that  the  rubrene  crystals  as  grown  are  epitaxially 
locked  in  the  crystallographic  direction  of  the  3 -fold  symmetric 
h-BN  crystal.  Furthermore,  from  the  comparison  between  the 
measured  GIXD  patterns  and  the  computed  diffraction  patterns 
for  a  series  of  0  angles,  we  extract  the  in-plane  crystallographic 
relationship  between  a-axes  of  rubrene  and  h-BN  with  an  offset 
of  ca.  3°  (see  Supporting  Information,  Figure  S2,  S3),  which  is 
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Figure  2.  VdW  epitaxial  relation  between  rubrene  and  h-BN  crystals,  a)  Bright-field  TEM  image  and  selected-area  electron  diffraction  pattern  of  the 
rubrene  crystal  grown  on  h-BN.  The  rubrene  (200),  (020)  and  h-BN  jloTo}  planes  are  indicated  by  orange,  yellow  and  green  circles,  respectively, 
b)  Schematic  of  GIXD  measurements,  c)  Measured  and  calculated  2D  GIXD  patterns  of  rubrene  crystals  grown  on  h-BN.  Green  and  white  circles 
represent  the  calculated  diffraction  patterns  for  orthorhombic  rubrene  and  h-BN  crystals,  respectively.  The  two  dominant  bright  features  observed  at 
the  diffraction  point,  (qn  qz)  ~  (2.4  A-1,  2.5  A-1)  and  ~  (1.2  A-1,  1.7  A-1)  originate  from  the  Si  substrate  underneath  the  thermally  grown  Si02  layer, 
d)  Representative  azimuthal  (b)-scan  of  GIXD  for  the  (256)  plane  of  rubrene,  which  shows  the  regular  diffraction  peaks  with  6-fold  symmetry,  e)  Top- 
view  image  of  rubrene  molecules  on  the  h-BN  surface  using  a  ball-stick  model.  The  black  box  defines  the  supercell  used  in  the  calculations  with  the 
definition  of  the  in-plane  misorientation  angle  (0).  f)  Plot  of  the  binding  energy  per  rubrene  molecule  as  a  function  of  0  obtained  from  first-principles 
calculations.  The  red  and  black  curves  represent  the  calculations  with  and  without  vdW  interactions,  respectively. 


similar  to  the  ca.  4°  offset  observed  in  the  SAED.  It  is  inter¬ 
esting  to  note  that  despite  the  relatively  weak  molecule-sub¬ 
strate  interactions  that  allow  the  structural  parameters  of  the 
rubrene  crystals  on  h-BN  to  be  close  to  those  of  free-standing 
rubrene  crystals,  the  interaction  is  strong  enough  to  provide 
epitaxial  registry  between  the  two  atomic  lattices. 

To  further  understand  this  epitaxial  relation,  we  per¬ 
formed  first-principles  binding-energy  calculations  using 
density  functional  theory  (see  Experimental  Section  and  Sup¬ 
porting  Information  for  more  details).  We  considered  super¬ 
cells  composed  of  a  single  layer  of  the  rubrene  crystal  on  top 
of  a  h-BN  layer  at  different  in-plane  misorientation  angles  (0) 
between  the  a-axes  of  rubrene  and  h-BN.  Figure  2f  shows  the 
total  energy  per  rubrene  molecule  as  a  function  of  0  using 
two  different  density  functionals  with  and  without  vdW  inter¬ 
actions.  In  both  cases,  the  minimum  energy  configuration  is 
obtained  when  small  misorientation  angles  of  less  than  10°  are 
present.  In  particular,  when  vdW  interactions  are  taken  into 
account  between  rubrene  and  h-BN  layers,  the  misorientation 
angle  becomes  even  smaller  to  6°  and  the  system  is  more  sta¬ 
bilized  with  reduction  of  energy  by  0.20  eV  per  rubrene  mol¬ 
ecule.  This  agreement  with  experimental  observations  strongly 
suggests  that  a  weak  vdW  interaction  can  still  lead  to  epitaxial 
growth  of  rubrene  crystals  on  h-BN  substrates. 

The  vdW  epitaxy  of  rubrene  films  allows  us  to  achieve  large 
single-crystalline  domains  on  single-crystal  h-BN.  Figure  3 


displays  a  polarized  optical  microscopy  image  for  a  typical 
rubrene  crystal  grown  on  h-BN.  The  contrast  in  the  polarized 
images  at  different  angles  is  dependent  on  the  in-plane  crys¬ 
tallographic  orientations  due  to  the  anisotropy  of  the  rubrene 
crystal. [17]  The  domain  size  is  determined  by  the  size  of  the 
h-BN  flakes,  indicating  that  the  rubrene  is  a  single  domain. 
The  large  single-domain  growth  of  epitaxial  films  is  critical 
to  achieve  intrinsic  device  characteristics  with  high  carrier 
mobility,  since  grain  boundaries  can  significantly  decrease  the 
mobility  and  thus  deteriorate  overall  device  performance  s.[26] 

Figure  4  displays  how  we  create  field-effect  device  struc¬ 
tures  on  these  epitaxially  grown  rubrene  crystals.  Electrodes 
are  required  not  only  to  contact  the  channel-dielectric  interface 
between  rubrene  and  h-BN  where  the  hole  accumulation  layer 
is  induced  in  our  silicon  back-gate  devices,  but  also  to  provide 
a  less  disruptive  surface  for  the  crystal  growth.  We  employ  vdW 
graphene  electrodes  that  serve  as  an  electrical  contact^27’28!  as 
well  as  an  epitaxial  growth  template  J23]  Figure  4a  shows  a  sche¬ 
matic  diagram  of  the  procedures  of  device  fabrication.  Prior  to 
rubrene  growth,  chemical-vapor-deposited  graphene  films  were 
transferred  onto  the  top  of  h-BN  layers  using  the  pre-patterned 
poly-dimethyl  siloxane  (PDMS)  stamp[27]  to  create  a  lateral 
graphene/h-BN/graphene  heterostructure.  During  this  transfer 
process,  the  h-BN  surface  over  the  channel  region  does  not 
contact  the  PDMS  and  therefore  remains  clean,  as  confirmed 
by  AFM  imaging  (Figure  4b).  The  separation  between  the  two 
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Figure  3.  Rubrene  crystalline  films  with  large  single-crystalline  domains  grown  on  h-BN.  Polarization  optical  characterization  for  two  representative 
samples.  The  normal  (1st  column)  and  polarized  (2nd-4th  columns)  images  at  different  angles  are  shown  for  each  samples.  The  last  column  shows 
the  grain  map  overlapped  on  the  optical  image  of  h-BN  flakes.  Crystal  grains  and  grain  boundaries  are  distinguished  clearly  by  an  optical  contrast  in 
polarized  images  and  indicated  by  dashed  lines  with  different  colors  in  maps. 

graphene  electrodes,  corresponding  to  the  channel  length,  was  Atomically  thin  graphene  provides  an  optimal  2D  mor- 
10  pm,  and  channel  width  determined  by  the  width  of  the  h-BN  phology  for  the  rubrene  crystal  growth  on  h-BN  and  is  less  dis- 
flakes  was  typically  in  the  range  of  10-30  pm.  ruptive  to  the  surface  than  typical  three-dimensional  (3D)  metal 
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Figure  4.  Fabrication  and  electrical  characterization  of  FETs  using  graphene  electrodes,  a)  Schematic  illustration  for  the  fabrication  of  graphene  elec¬ 
trodes  on  h-BN.  b)  AFM  image  of  the  patterned  CVD-graphene  electrodes  on  the  h-BN/Si02  substrate,  c)  Optical  images  of  the  device  before  (top) 
and  after  (bottom)  rubrene  crystal  growth  on  the  stacked  graphene/h-BN  structure.  d,e)  Representative  output  (d)  and  transfer  (e)  characteristic 
curves  of  FETs  with  graphene  electrodes,  f)  Histogram  of  mobility  measured  from  27  devices.  The  average  and  maximum  values  of  mobility  are  5.1  ± 
2.7  cm2  V-1  s_1  and  11.5  cm2  V-1  s-1,  respectively. 
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electrodes.  For  example,  when  Au  metal  electrodes  were  used, 
small-domain  polycrystalline  films  were  grown  on  the  contact 
area  presumably  due  to  random  heterogeneous  nucleation 
along  the  3D  edges  as  wells  as  lack  of  preferential  assembly  of 
molecules.  This  results  in  structural  disorders  such  as  grain 
boundaries  that  hindered  carrier  injection  into  the  channel, 
causing  extremely  non-linear  current-voltage  curves  with 
large  contact  resistances  of  a  few  hundred  mega-ohms  (see 
Supporting  Information,  Figure  S4).  In  contrast,  as  shown  in 
Figure  4c,  the  high-quality  rubrene  films  were  grown  uniformly 
across  the  graphene/h-BN/graphene  hetero structure  where 
graphene  only  has  an  one-dimensional  edge  with  the  rubrene 
crystals.  We  also  find  that  output  characteristics  of  field-effect 
transistors  (FETs)  are  significantly  improved  compared  with 
those  of  the  devices  with  metal  electrodes.  In  particular,  the 
series  resistance  at  low  source-drain  bias  (<1  V)  was  signifi¬ 
cantly  reduced  down  to  a  few  MQ  (Figure  4d),  an  improvement 
of  two  orders  of  magnitude  when  compared  with  Au  electrodes. 
Note  that  the  contact  resistances  (about  several  kQ  cm)  are  still 
high  compared  with  those  of  single-crystal  devices^29!  because 
our  epitaxial  films  can  be  structurally  disordered  in  the  vicinity 
of  randomly  stacked  h-BN/graphene  boundaries  (see  Sup¬ 
porting  Information,  Figure  S5,  S6). 

The  aforementioned,  graphene-contacted,  vdW  heterostruc- 
tured  devices  exhibited  excellent  device  performances  as  organic 
transistors  despite  an  imperfect  contact.  Figure  4e  shows  the 
transfer  characteristic  curves  of  FETs  with  negligibly  small  hys¬ 
teresis  and  high  on/off  ratio  of  ca.  106,  indicating  very  low  den¬ 
sity  of  deep-level  trap  states  at  the  interface.  In  addition,  field- 
effect  mobility  in  the  saturation  regime^30!  was  substantially  high; 
the  average  and  maximum  values  were  5.1  ±  2.7  cm2  V-1  s-1  and 
11.5  cm2  V-1  s_1,  respectively.  Figure  4f  shows  a  histogram  of 
mobility  measured  from  27  devices.  Note  that  the  linear  mobility 
was  typically  an  order  of  magnitude  lower  (0.1-1  cm2  V-1  s-1) 
than  the  mobility  we  extracted  from  the  saturation  regime  (see 
Supporting  Information,  Figure  S6).  Nevertheless,  the  highest 
mobility  observed  in  these  epitaxial  vdW  heterostructures  is 
comparable  to  those  of  free-standing  single  crystals  and  repre¬ 
sents  a  record  for  organic  films  grown  on  any  substrates  (see 
Supporting  Information,  Table  SI).  It  should  be  also  noted  that 
the  substantial  variation  in  mobility  could  be  attributed  to  the 
anisotropy  of  mobility  according  to  epitaxial  crystal  orientations 
determined  randomly  by  exfoliated  h-BN  crystals,!15!  rather  than 
grain  boundaries  that  typically  degrade  carrier  mobility  by  more 
than  an  order  of  magnitude.!26!  All  these  device  characteristics 
including  negligible  hysteresis  and  high  mobility  strongly  sup¬ 
port  the  formation  of  highly  ordered  crystalline  films  with  an 
electrically  clean  semiconductor/dielectric  interface  that  is  con¬ 
sistent  with  our  structural  analyses. 

In  summary,  we  demonstrate  that  epitaxial  vdW  heterostruc¬ 
tures  of  organic  and  vdW  layered  materials  create  high-perfor¬ 
mance  organic  electronic  devices.  High-quality  rubrene  films 
with  large  single-crystalline  domains  were  grown  on  h-BN  dielec¬ 
tric  layers  via  vdW  epitaxy.  High  carrier  mobility  comparable  to 
free-standing  single-crystal  counterparts  was  achieved  by  forming 
interfacial  electrical  contacts  with  graphene  electrodes.  This  high¬ 
lights  the  advantages  of  heterostructured  vdW  dielectrics  and 
contacts  over  commonly  used  materials  to  achieve  high-perfor¬ 
mance  organic  electronic  devices.  More  generally,  the  epitaxial 


approach  to  fabricate  organic/inorganic  vdW  heterostructures 
can  be  readily  expanded  to  numerous  other  organic  and  layered 
materials  for  various  electronic  and  optoelectronic  applications. 


Experimental  Section 

Material  Growth:  Rubrene  crystalline  films  were  grown  by  a  vapor- 
phase  transport  method.  A  source  of  rubrene  powder  (ca.  1  mg),  which 
was  purified  two  times,  was  placed  in  the  sublimation  zone  at  230  °C. 
The  Si02/Si  substrate,  on  which  h-BN  flakes  were  exfoliated,  was  placed 
in  the  crystallization  zone  at  180  °C.  Prior  to  the  growth,  the  h-BN 
surface  was  cleaned  via  thermal  annealing  at  400  °C  under  vacuum  (ca. 
1  mTorr).  AFM  surface  characterization  after  annealing  reveals  that  we 
still  have  an  atomically  flat  and  clean  substrate  without  much  of  damage 
or  contamination.  During  the  growth,  the  chamber  was  under  vacuum 
(ca.  1  Torr),  with  flowing  the  high-purity  argon  gas  (99.999%)  as  a 
carrier  gas  with  the  flow  rate  of  200  seem.  The  growth  was  performed 
for  5  minutes,  and  then  the  system  was  cooled  down  slowly  to  room 
temperature. 

Material  Characterization :  For  the  TEM  measurements,  rubrene  films 
were  grown  directly  on  h-BN  flakes  suspended  onto  a  copper  TEM 
grid  using  a  mechanical  exfoliation  method.  Scanning-TEM  images 
and  SAED  patterns  were  obtained  at  a  low  operation  voltage  of  80  kV 
to  avoid  the  electron-beam  damage  with  the  TEM  (JEOL,  JEM-100CX). 
For  the  GIXD  measurements,  a  single  piece  of  an  h-BN  crystallite  with 
the  lateral  dimensions  larger  than  100  pm  was  transferred  onto  a  Si02/ 
Si  substrate  using  the  PDMS  stamping  technique,  and  the  rubrene 
crystals  were  grown  the  following  the  method  described  above.  This 
allows  us  to  obtain  GIXD  patterns  from  the  rubrene  films  grown  on 
a  single  crystallite  of  h-BN.  The  experiments  were  performed  at  the 
National  Synchrotron  Light  Source  (NSLS)  of  Brookhaven  National 
Laboratory  on  the  undulator-based  X9  beam  line,  X-ray  photons  of 
14.0  keV  (0.08856  nm).  The  incident  beam  has  a  grazing  incidence  angle 
with  the  sample  surface.  The  X-ray  beam  was  focused  to  100  pm  in  the 
horizontal  and  50  pm  in  the  vertical  direction.  The  beam  was  aligned 
with  the  sample  using  a  phosphor  to  calibrate  beam  position;  the 
sample  was  laser-levelled  to  the  in-plane  rotation  axis  and  rotated  in  15° 
steps  to  determine  in-plane  registry  with  the  underlying  h-BN  substrate. 
Additional  details  are  provided  in  the  Supporting  Information. 

First-Principles  Calculations:  The  simulations  reported  here  are  based 
on  density-functional-theory  calculations  using  the  SIESTA  code.!31! 
The  generalized  gradient  approximation!32!  and  nonlocal  van  der  Waals 
density  functional!33!  was  used  together  with  double-q  plus  polarized 
basis  set,  norm-conserving  Troullier-Martins  pseudopotentials!34!  and 
a  mesh  cutoff  of  150  Ry.  We  used  a  Fermi-Dirac  distribution  with  an 
electronic  temperature  of  kBT  =  21  meV.  Additional  details  are  provided 
in  the  Supporting  Information. 


Supporting  Information 

Supporting  Information  is  available  from  the  Wiley  Online  Library  or 
from  the  author. 
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